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Abstract 

In the inertial range of fully developed turbulence, we model the vertex 
network dynamics by an iterated unimodular map having the universal 
behavior. Inertial range anomalous scaling for the pair correlation func- 
tions of the velocity and the local energy dissipation is established as a 
consequence of the chaotic behavior of the unimodular map when the 
Feigenbaum attractor looses stability. The anomalous corrections deter- 
mined by the Feigenbaum constant 77 to the Kolmogorov's spectra are 
larger than those observed in experiments. 

Turbulence Modeling and Turbulence Control in fluid dynamics 
is a paradigmatic problem and successful techniques in this domain 
may serve as models to deal with other complex extended systems. 
It is generally accepted now that a real understanding of turbulence 
requires calculations and models that are globally sensitive to all 
length scales. 

To elucidate this conviction, let us briefly explain the present 
situation in the quantum fleld theory of fully developed turbulence 
of incompressible fluid which has been developed for the last two 
decades (see for a review). In order to be speciflc we consider the 
correlator K = (f (k)f (— k)) of the velocity fleld v in the momentum 
representation {k = |k| and (. . .) denote the time averaging). It is 
found from the Dyson equation 

K-^ = uk'^ + -j:{k), (1) 
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where m is the inverse integral turbulent scale, u is the viscosity 
parameter, and is an infinite sum of all 1-irreducible diagrams 
of the Wyld's diagrammatic technique [0]. In the framework of 4— 2e 
expansion the series (|I|) takes the form |^ 



2 + m2) 



1 + Y,(gk-^Tcn{m/k, 



n=l 



(2) 



The inverse viscous dissipative length A plays the role of a maximum 
momentum in the problem. The Reynolds number Re= {X/m}^^^ is 
assumed to be large. For k ^ m <^ X and e > 0, the dimensionless 
parameter of the expansion gk'"^^ ~ [k/X)^"^^ /v in is not small 
and it is necessary to sum the series. This problem reduces to a de- 
termination of the asymptotic value of = K{xk, x^^) fo^^ X 
(everything is fixed except for x) . It is solved by the renormalization 
group procedure (KG) for any finite ratio m/k 

However, the coefficient c„ in (|^) has additional singularities as 
m/k which cannot be handled by RG since they have nothing 
to do with the scale invariance property (the parameter m does 
not involve the renormalization). These singularities express the 
intermittency phenomenon (0, Chap. 25) that is the actual reason 
why the bare RG is ineffective in turbulence [0]. 

The m/k asymptotics of the series (0) is found due to the 
Short Distance Expansion method 0, introduced by Wilson. It 
leads to the form 



K{k) = Ak 



-11/3 



(3) 



where A is the Kolmogorov's constant, 6j is the coefficient analytic 
in m/k, and Ap. is the critical exponent of a set of composite oper- 
ators Fi. The series in (Q) would give the corrections to the leading 
Kolmogorov's asymptotics if all Ap^ > 0. 

A number of works has been devoted to the computation of Ap^ 
for various sets of composite operators in hydrodynamics , and 
magnetohydrodynamics |T^. It is found that there are infinitely 
many composite operators for which Ap^ < 0. An example of such 
a "dangerous" one is the composite operator of local energy dissi- 
pation E = —u{diVk + dkViY /2 with A^; = 0. The contributions 
to (Rf) due to such operators are not small and it is necessary to 
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sum the infinite series in (^). There is no rehable technique for 
such a summation found. To proceed further, remaining of the idea 
of homogeneously fully developed turbulence, one must assign the 
statistics to each set of composite operators with Ap. < 0. 

It is conventional now that one has to employ a model, in which 
all scales are presented equally to achieve the further progress. In 
this Letter we suggest such a model. 

The hydrodynamics of the incompressible fluid is described by 
the Navier-Stokes equation 

VtV^ = uAvi, Vt = dt + {vd), (4) 

where fj(x, t) is the transverse {diVi = 0) velocity field, z/ is the 
viscosity. To meet the energy balance in the system, one has to 
compensate the viscous dissipation by an energy pump of power W 
with some spectral density ^{k), 

W = — ^ [ dk d{k). (5) 



(2vr)3 , 

In and in the descendant papers 0-|lO[, the pump source 

has been introduced in (^) by the Gaussian white noise /. The 
power model for the spectral density d{k) ~ k{k'^ +m^)~^ chosen in 
literature defines the Gaussian statistics for / completely, {f f){k) ~ 
Nevertheless, such a model is not uniquely determined: One 
can introduce an arbitrary bounded value function \){rn/k) (f)(0) = 
1) into the Gaussian covariance. This function did not involve the 
renormalization procedure in From the physical point of 

view, this fact means that the field / has some redundant degrees 
of freedom which do not change to the scaling ones in the process 
of renormalization. 

In actual experiments, the power W pumped into the inertial 
range of scales m -C <^ A is delivered by the large scale vertices 
entering the system sometime in a variety of sizes {k ~ m), energies, 
and enstrophies. These vertices dissolve due to the strong nonlin- 
earity of hydrodynamical interactions into ones of the smaller scales 
{k > m) forming the direct energy cascade. Alternatively, a small 
scale vertex would associate with others in the fractal sets contribut- 
ing towards the inverse cascade. The resulting picture constitutes a 
vortex network encoded, at each moment of time, by a discrete set 
of ratios < {m/fc}„ < 1. This network has much in common with 
dynamical systems. 
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We shall consider the pumping generated by some dynamical 
system, 

[t/r] 
n=0 

Here r > is a time step, the square brackets [. . .] denote the 
integer part, the iterates of a 1— parameter unimodular map 

ni T : X ^ X of the phase space X = [0, 1]. We also suppose that 



the map T possesses the invariant measure /xt- 

As an example of such a map, one can consider the triangular map 
Tz = r(l — 2|i — or the logistic map Tz = rz{l — z) which is 
known to describe the angles of a strongly damped kicked rotator. 
However, the existence of invariant measure for the latter map is 
not proven for any r except r = 4. The dynamics of unimodular 
maps has been studied extensively last decades. For the values of 
the control parameter 1 < r < Too, the relevant Liapunov exponent 
is always negative (except the bifurcation points r„ when it becomes 
zero), while for < r this exponent is mostly positive indicating 
chaotic behavior. 

The chaotic behavior observed in the unimodular maps results 
from the successive pitchfork bifurcations, which provides the mech- 
anism for the successive doubling of the fixed points. It is important 
that the correlations between the j-shifted sequences 



C{j) = lini 





n>Q 



[t-t'/r] 



^nZn+j (7) 



_t-t'_ 

generated by the unimodular map decay with a power law in j as 



0, 

C{j)c^r''C{l,p) + 0{r-r^). (8) 

Here t] is the universal constant ||13[, p is the "scaling function" 
invariant with respect to the doubling transformation T^. 

An external force fr{'x,t) : W^^^ x X ^ M'^"*"^, which we intro- 
duce in d^) is regarded as a priori arbitrary chaotic process Wr = 
T"^/^ J2n>o ^"^(^ ~ ^^)^ G X!, generated by a deterministic chaotic 
dynamics of the map T {d is the space dimensionality). 

In the language of quantum field theory, the time step r > plays 
the role of the limiting time scale, dividing the statistics of "fast" 
{t < t) and "slow" {t > r) modes. It is evident that the statistics 
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of the fast modes f (x, t < r) is not affected by the dynamics of the 
map T, so that one can treat as a Gaussian white noise with the 
covariance 



A.s(x-y,t-t') = (27r)"^ dio dk Prs^{k)e 



,i[k(x-y)-c^{t-t')] 



(9) 



where Psr = Ssr — kgkr/k'^ is the transversal projector. Following 
we chose the spectral density in the power form d{k) = gv^k^~'^~'^^ , 
in which is a coupling constant. The actual value of the parameter 
of regular expansion is e = 2, so that D(/c) represents a power model 
for 5(k) located in the infrared {k ~ 0) region. The other way 
round, the statistics of the slow modes fr(x, t > r) is essentially 
non-Gaussian. 

Then the velocity field correlator K = (v(x, t)v(y, t')) factors in 



K 



Ki^-y,t-t')J2Crn{t-t';qi,...qn), t > t' . (10) 



n>2 



Here .^(x — y, t — t') is the correlator in the fast modes Gaussian 
statistics, and the infinite sum describes the chaotic process as t > r. 
All correlations 



C 



lim 



t-t' 



[t-t'/r] 
p>0 



^P+qn! Ifc — 



t — tk 



e z. 



(11) 

{t > . . . > tk > tfc-i > . . . t' > 0) exist since there is the continuous 
invariant measure fiT- 

The correlator (|T^) is finite for any ratio m/k fixed: K is found 
from the Dyson equation (|l]) in the framework of some quantum field 
theory with the multiplicatively renormalized action functional 



^1 



R 



2 




dk dt \'D\' 




+ dkdt v'[-dtv - {vd)\r + z/Z^Av] 



(12) 



of the basic v and auxiliary v' fields. D is the covariance (^. The 
only renormalization constant z/q = uZiy is required to subtract sin- 
gularities, dJl is the renormalization mass parameter. All correlation 
functions of the renormalized quantum field theory are finite for any 
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fixed m/k. In the large-scale asymptotics s = /c/9Jl— s>0, in(i = 3, 
they demonstrate the scaling behavior with the Kolmogorov's criti- 
cal indices of velocity = —1/3 and time = —2/3. As we have 
mentioned above, in the limit m/k 0, the additional singularities 
arise in ^. We demonstrate, however, that (|IU]) is still finite. 

In the majority of experimental data, the correlator K and the 
correlation function of the local energy dissipation E = ($(x)$(y)), 
$(x) = 5jfj(x)9jfj(x), have a power law behavior 

ir(A;) ~ A;-("/3+5)^ E(x, y) ~ |x - y]"^ (13) 
with small positive values of indices A ~ 0.2, 0.02 < 5 < 0.07 in 



accordance with [|J]'I1I11- The Kolmogorov's theory predicts A = 
5 = 0. Now we demonstrate that in the scaling limit r ^ 0, at 
Too (i.e., when the Feigenbaum attractor of the map T looses its 
stability) the infinite sum in ([l0|) provides the anomalous scaling 
A >(5 > 0. 

First, we show that the pair correlation provides the leading 
contribution to the sum in ([T0| ) oc (r/t — if:')''. The contributions 
due to the triple and quadruple correlations etc. decay much more 
rapidly as r — 0. To prove it in an elegant way, we use the theorem 
followed which is analogous to the Wick's theorem of quantum field 
theory: The average of the product is obtained as a sum of all 
possible paired averages. 

Theorem 

In the scaling region r ^ 0, the following expansion has place 

n 

C^„(t; qi, . . . q^) = ^ C2r(t; qi) W C2r(t; q^ - q^)- 

The elementary proof of the theorem is given in the Appendix. This 
theorem gives a key to diagram technique in dynamical systems. 
The diagrammatic approach to computation of arbitrary correlation 
functions in dynamical system theory has been developed recently 
One can estimate the scaling asymptotics for correlation 



m 



functions as Cm oc r'^'"'^ {n is assumed to be an even number) then 
the result on the scaling behavior seems obvious. 

Interesting in the scaling asymptotics r — > 0, we neglect all higher 
order correlations in (0) except C2T, so that, for r — > Too, in the 
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inertial range, in d = 3, (|T0|) reads as follows 

K{k, t-t') = A'k-''/' {t^)' + ^^"^ ~ ^^^^ 
in which A' is some modified Kolmogorov's constant. The asymp- 



totics (JT^ can be treated as a consequence of modification of the 



critical regime in the theory (|T^) as the result of rescaling: 

v = v^, v' = v'r'', ~9 = 9e^^ k = k^, i=tr', (is) 

in which 



2 + d-2A^' 

is the Kolmogorov's dimension of velocity. The critical exponents 
of all "-quantities are known (they equal to the Kolmogorov's values). 
The scaling parameter of the doubling transformation r does not 
involve the theory (|T^), so that it has no definite critical dimension. 



Therefore, the modification of scaling asymptotics s = k/DJl ^ 
due to (|l^) goes from the rescaling with respect to the dimensional 
parameter t cx s~^/^. The upper critical dimension (above which the 
hydrodynamical interaction becomes infrared unimportant, and the 
field theory is asymptotically free) for the rescaled model is changed 
to 4 — 4/3(2 + (i)(T and is not equal to the logarithmic dimension 
di = 4. In dimension d = 4 — 2e, the hydrodynamical interaction 
is important for 2/3(2 + ci) a < e (this condition is satisfied for the 
actual value e = 2). Thus, the quantities 

A'fc = 1 - 2/3a, a; = A, - 2/3a, A[ = At + 4/3a, (16) 

generated by (|I3]) are the entire critical dimensions of momentum, 
velocity, and time in turbulence launched by the dynamical system. 
The simple power counting in (|l^) gives the asymptotics 



K ~ /t~(ii/3+5j^ ^ 2ri/3a ~ 0.2014722 .... 

The renormalization of the local energy dissipation operator E with 
mixing has been given in [^] in details. It was shown that the 
marginal contribution into asymptotics is risen due to the linear 
combinations of composite operators including (vAv)(x), which has 
the critical exponent A^; = 4 — S'jt = 0, where 7t = 2 + At = 4/3 
is the anomalous dimension of time in the Kolmogorov's theory. In 
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the rescaled theory, the anomalous dimension of time is changed to 
7^ = 2 + ~ 0.7107451 . . . , so that the resulting critical dimension 
~ 1.8677646. . . > is no more marginal. This corresponds to 
the power law behavior 

E(x,y) ~ |x-y|"^, A ~ 1.1322354.... 

In conclusion, we have studied the model of turbulence driven by 
the deterministic chaotic dynamics generated by a unimodular iter- 
ated map. The critical scaling is established when the Feigenbaum 
attractor looses stability. When the new time scale comes into play, 
the Kolmogorov's critical regime is modified and the anomalous ex- 
ponents becomes positive X > S > 0. However, their absolute values 
prescribed by the Feigenbaum constant rj are larger than those in 

mm- 



Appendix 

A proof sketch of the Theorem. 

Let us consider the product of partial sums 

^^(q, qi • • • qn) 

= ~ ZpZp^q— ^Pi + qi^pi+q2 ' ' ' ~ ^Pn + qn-i^Pn + qn- 

P Pi pn 

Rename the counting variable in each sum, for example 



ra+qi 

2^r^r+{q2-qi)- 



1 \ - _ n + qi 1 X - 

~ / ^ ^Pi + qi^pi+q2 ~ ~~~ / , 

pi ' r=l+qi 

As 71 — > oo, this partial sum tends to C2(q2 — qi)- It remains to note 
that 

Cn(qi,---qn) = ^ 5ppi---5pp„_i*P(q,qi...qn). 

qi...qn-i 
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